Majority-vote model on hyperbolic lattices 
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We study the critical properties of a non-equilibrium statistical model, the majority-vote model, 
on heptagonal and dual heptagonal lattices. Such lattices have the special feature that they only can 
be embedded in negatively curved surfaces. We find, by using Monte Carlo simulations and finite-size 
analysis, that the critical exponents 1/u, j3/v and j/u are different from those of the majority- vote 
model on regular lattices with periodic boundary condition, which belongs to the same universality 
class as the equilibrium fsing model. The exponents are also from those of the fsing model on a 
hyperbolic lattice. We argue that the disagreement is caused by the effective dimensionality of the 
hyperbolic lattices. By comparative studies, we find that the critical exponents of the majority-vote 
model on hyperbolic lattices satisfy the hyperscaling relation 2/3 /v + y/u — D a g, where D c s is an 
effective dimension of the lattice. We also investigate the effect of boundary nodes on the ordering 
process of the model. 

PACS numbers: 05.50.+q, 02.40.Ky, 64.60.Cn, 05.70.Ln 
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I. INTRODUCTION 

Recently, there has been a growing interest the crit- 
ical behavior of statistical-physics models on curved 
surfaces — ranging from spin models, such as the ferro- 
magnetic Ising model J3, H, H, 0] , the XY model 
the Hcisenberg model Q , the g-state clock models 
to other traditional models, such as percolation 
fusion [11], etc. One reason for this interest is that many 
newly discovered soft materials (e.g., carbon nanotubes) 
show a negatively curved structure in the nanoscale [l^ | . 
One peculiar feature of negatively curved surfaces is that 
their boundary is a finite fraction of the whole geome- 
try This structure has been verified having a non- 
trivial impact on the critical behavior of many models 
of statistical physics. For example, in the context of 
the Ising model, significant shifts in static and dynamic 
critical exponents toward the mean-field values were no- 
ticed [H, Q ; small-sized ferromagnetic domains were ob- 
served to exist at temperatures far greater than the crit- 
ical temperature An apparent zero-temperature ori- 
cntational glass transition in the XY spin model on a 
negatively curved surface was recently demonstrated Q. 

These findings motivate us to investigate another spin 
models on neg atively curved surfaces. The majority- vote 
model [II [ll [II El is a simple non-equilibrium model 
exhibiting up-down symmetry that has been argued to 
belong to the universality class of the equilibrium Ising 
model 18]. Oliveira first verified this conjecture on a 
square lattice with periodic boundary conditions (i.e., 
a torus) (l4| . Subsequently, the majority- vote model 
has been investigated on regular lattices (with dimension 
larger than two) [H, H3, H3| , random lattice [HI , di- 
rected or undirected random graphs 0, HH , small world 



networks [24|, and scale- free networks [25(, etc. Very re- 
cently, it has been found that the critical behavior of the 
majority-vote model on square lattice is also indepen- 
dent of transition rates (e.g., the Glauber or Metropolis 
rates) [19]. It has also been observed that the majority- 
vote models defined on different complex networks belong 
to different universality classes [22|, 0, [24|, [25| . 

Our goal in this contribution is therefore to identify 
the critical behavior of the majority-vote model when 
the underlying lattice is embedded in a hyperbolic sur- 
face, in particular the heptagonal and the dual heptago- 
nal lattices, and investigate if such an interaction struc- 
ture is capable of modifying the critical exponents. To 
this end, we use Monte Carlo (MC) simulations and stan- 
dard finite-size scaling techniques to determine the crit- 
ical noise parameter q c (the main control parameter of 
the majority-vote model, as well as the critical exponents 
1/v, fl/v and yjv. Our numerical results suggest that 
the critical exponents, in the stationary state, are differ- 
ent from those of the Ising model confined to regular and 
hyperbolic lattices. 

In the following Sec. [Til we define our model, describe 
the quantities we measure and the computational details. 
In Sec. IIIH we present our numerical results and analysis. 
Finally, we summarize and contextualize the observations 
in Sec. [TV] 
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II. MODEL AND SIMULATION 
A. Majority-vote model 

Following Refs. 0, EE EE E3, EH, [H, H Hi, the 

majority-vote model with noise is defined by a set of spin 
variables where each spin is associated to one node 
of the heptagonal lattice and can take the values ±1. 
The system evolves as follows: For each spin i, we first 
determine the majority spin of i's neighborhood. With 
probability q the i takes the opposite sign of the majority 
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spin, otherwise it takes the same spin as the majority 
spin. The probability q is called the noise parameter and 
plays the same role of temperature in equilibrium spin 
systems. In terms of q, the probability of a single-spin- 
flip is given by 



the Binder's fourth-order cumulant U. These quantities 
arc defined as follows: 
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where S(x) — sgn(a;) if x ^ and 5(0) = 0, and the sum- 
mation is over all the neighboring spins of the focal site i. 
The transition probability fl} satisfies up-down symme- 
try under the change of signs of the spins in the neighbor- 
hood of i. In the limit of zero noise, the majority-vote 
model is identical to the ferromagnetic Ising model at 
zero temperature [HI, [H| . 



B. Heptagonal lattice and dual heptagonal lattice 

Figure [1] shows two examples of the heptagonal lat- 
tice and the dual heptagonal lattice. One peculiar prop- 
erty of the heptagonal lattice is that, if we consider 
the innermost heptagon as the level one, then the num- 
ber of nodes of a heptagonal lattice with level I can 
be calculated by using the formulation 0, 0, N(l) = 

1 + ELiK^T^P " i^T^yi wh ic h grows exponen- 



tially with level I. In other words, the ratio of the perime- 
ter to the area of the lattice remains finite (about 0.62) 
in the thermodynamic limit I — ► oo [ij . If we make map 
each heptagon in the heptagonal lattice to a node, and 
put a link between adjacent heptagons, then we get the 
dual heptagonal lattice [Fig. [1Kb)]. Since these lattices 
can only be embedded in a hyperbolic surface with a con- 
stant negative curvature, they are also called hyperbolic 
lattices. 





FIG. 1: A heptagonal lattice with level 3 (left), and the 
dual lattice with level 4 (right). For the heptagonal lattices 
with level I — 2,3, 4, 5, 6, 7, the total number of nodes are 
35, 112, 315, 847, 2240, 5887, respectively. For the dual lattice 
with level I — 3, 4, 5, 6, 7, 8, the total number of nodes are 
29, 85, 232, 617, 1625, 4264, respectively. 



C. Measurements and finite size scaling 

In order to study the critical behavior of the model we 
consider the magnetization M, the susceptibility Xi and 
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where N is the total number of nodes of the hyperbolic 
lattice, (...)t denotes thermodynamics averages taken in 
the stationary regime, and {■ ■ - )c stands for configura- 
tional averages. 

The above quantities are functions of q and N, in the 
critical region, we expect the following finite-size scaling 
relations 



M(N, q) = N-^' v M{N l l u e) 
X(N,q)=N^x(N 1 ^e) 
U\N,q) = N 1 I v U'{N 1 I"e) 



(5) 
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where e = q — q c and U' is the derivative of Binder's 
fourth-order cumulant with respect to the noise. By the 
standard finite-size scaling approach [27| , we assume scal- 
ing functions M, x, and U that are continuous and diffcr- 
entiable in the vicinity of the critical noise q c . From the 
size dependence of M and % we can obtain the exponents 
j3/u and j/p, respectively. One alternative way to detect 
j/v, since it also scales as N 1 ^ [27l ]. is by plotting the 
maximum value of the susceptibility versus TV. 



D. Computational method 

We implement our MC simulations on the hyperbolic 
lattices with various system sizes starting with all spins 
up and going from low noise to high noise. It has been 
pointed out that this method can reduce the relaxation 
time considerably, especially in the low noise limit, com- 
pared to starting with spins randomly oriented up or 
down [26] . For each given q, we simulated systems of size 
N = 112,315,847,2240,5887 for the heptagonal lattice, 
and N = 85, 232, 617, 1625, 4264 for the dual heptago- 
nal lattice. For the sake of comparison, we also studied 
the majority-vote on square lattice with free boundary 
condition. In all simulations, we first wait 90 000 MC 
steps to let the system attain stationary state (in the 
high noise case, this number reduces to 50 000), and then 
ran another 30 000 MC steps to get the average values. 
One MC step contains a sweep of the spins in a ran- 
dom sequence. After every tenth MC steps, we reshuffle 
the random sequence. Each data point presented below 
are averages over 500, 300 and 200 trials for N < 1000, 
1000 < N < 4000, and N > 4000, respectively. 
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FIG. 2: (Color online) Majority- vote model on a two dimensional square lattice with a free boundary condition. Magnetization 
M (a), susceptibility \ (b), and reduced fourth-order cumulant U (c), as a function of the noise parameter q for several values 
of the system size N. In (c), within the accuracy of the data, all curves intersect at q c = 0.074(5). 
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FIG. 3: (Color online) (a) Data collapse of the magnetization M shown in Fig. [2ja). The exponents used for the data collapse 
are l/v = 0.5 and (3/v = 0.0625. (b) Log-log plot of the maximum of the susceptibility as a function of N. From it we estimate 
the critical exponent 7/1/ = 0.824(7) as the best fit of the data points, (c) Data collapse of the susceptibility shown in Fig.[2jc). 
The fitting exponents are p/v = 0.0625 and j3/v = 0.825. 



III. NUMERICAL RESULTS AND FINITE SIZE 
SCALING ANALYSIS 

Previous investigations have shown that the majority- 
vote model undergoes a phase transition from an ordered 
to a disordered phase at a critical value of cu This criti- 
cal value depends on the lattice topology [3, [HI, . It 
is worth noting that almost all previous studies assume 
a periodic boundary condition of the underlying topol- 
ogy. There are no previous results for the majority-vote 
model on regular lattice with a free boundary condition 
which is the case for the heptagonal and dual heptagonal 
lattices (Fig. [I}. Therefore we proceed to investigate the 
majority-vote model on a square lattice with free bound- 
ary condition. 



A. Majority-vote model on square lattice with free 
boundary condition 

In Fig. [3] we present MC results of the majority-vote 
model on square lattice with free boundary condition 
and plot the magnetization, susceptibility, and reduced 
fourth-order cumulant as functions of the noise param- 
eter q for several values of N. As can be noticed in 



Fig. [21 a) there is a phase transition from an ordered state 
(M > 0) to a disordered state (M wO). In Fig. [2(b), the 
susceptibility \ reaches a maximum in the critical re- 
gion for different system sizes which is another typical 
signature for the onset of criticality. The critical point 
q c can be detected from Fig. [2(c), where the curves of 
the reduced forth-order cumulant U for different N in- 
tersect with each other. We obtain q c — 0.074(5), which 
agrees quite well with that for the majority-vote model 
on a square lattice with a periodic boundary condition 

q c « 0.075 m m. 

The results of Fig. &&c) indicate that a vanishing frac- 
tion of boundary nodes do not change the critical value 
of q. In light of this point, we may expect the critical 
exponents to be the same. To check this assumption, 
we plot M N$l v versus (q - q^N 1 ^ in Fig f^a) using 
the critical exponents \jv = 0.5, fl/v — 0.0625 of two 
dimensional Ising model. (Note that if we measure the 
critical exponents in terms of linear dimension L instead 
of N = L 2 , then 1/^=1, /3/u = 0.125.) The excellent 
collapse of the curves for six different system sizes corrob- 
orates the estimations for q c and the critical exponents 
1/v and and verifies the Ising universality class of 
the phase transition. On the other hand, the curves for 
yJSS" 1 ! 1 ' versus (q — q^N 1 /" do not overlap in the critical 
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region if j/v = 0.875 is used (results not shown), which 
hints that the free boundary condition induces a strong 
finite-size effect on the fluctuation of average magnetiza- 
tion (the critical exponent 7/V). For this reason, we use 
finite-size scaling by plotting Xmax (the maximum value 
of the finite-size susceptibility) as a function of N. We 
get the slope 7/1/ = 0.824(7) of the best-fit line, as dis- 
played in Fig. [3](b) . The collapse of the curves for the 
rescaled susceptibility verifies this estimate [Fig. 0K C )]- 



B. Majority- vote model on heptagonal lattice 

We now turn our attention to the majority-vote model 
on the heptagonal lattice. The simulation results for M, 
X, and U as functions of the noise parameter q for dif- 
ferent N are summarized in Fig |4j As was shown in 
Fig. UK a), above the critical noise level, the magnetization 
disappears for larger system sizes, whereas it reaches a 
finite value in the subcritical region [Fig. HJa)]. The sus- 
ceptibility values in the critical region reach their maxi- 
mal values [Fig. 0Kb)], and the reduced fourth-order cu- 
mulants cross at the critical point, giving q c = 0.034(8). 

In order to study the universality class of the model, 
we proceed to measure the critical exponents. From 
Ref. [29(, we know that the critical exponent v can be 
obtained from the scaling behavior of certain thermo- 
dynamic derivatives with respect to q (for example, the 
derivative of the cumulant and the logarithmic deriva- 
tives of M). In Fig. 0Ka), we plot the maximum value 
of these derivatives as functions of system size on a log- 
log scale. 1/v obtained from these fits can be seen in 
Fig-EKa). From combining these two estimates, we obtain 
1/v = 0.30(2). In Fig. [5Kb), we display the JV-dependence 
of the magnetization at q c . From the slope of the dashed 
line, which corresponds to the best fit to the data points, 
we estimate the corresponding value of the critical expo- 
nent to [3/v = 0.114(5). Using these values we proceed 
to plot UN 13 '" against {q - q^N 1 ^. Fr om the finite- 
size scaling assumption that M is a universal function, 
we should, for the correct choices of 1/v and /3/v, find 
a data collapse in the critical region [28|. In Fig. 0Kc) 
we show that such a data collapse does indeed occur. 
In a similar way, we can determine the value of j/v by 
fitting the data for Xmax as a function of N in a log-log 
scale, whose slope predicts j/v = 0.721(7) [Fig.0Kd)]. By 
plotting xN 1 / 1 ' versus (q — q c )N 1 / V with 1/v = 0.3 and 
7/V = 0.722, however, we only get good collapse for the 
curves in the supercritical region, i.e., q > q c [Fig. 0[e)]. 
In the subcritical region the curves deviate, suggesting an 
anomalous scaling behavior. This effect can also be seen 
in Fig. 0Jd) for small system sizes. Since the boundary 
vanishes with size in Fig. 01 as does the deviation from 
the scaling collapse, but the boundary does not vanish in 
Fig. 0] and neither does the deviation, we conclude that 
the boundary is probably causing the anomalous scaling 
behavior. 



C. Majority-vote model on dual heptagonal lattice 

We proceed to investigating the majority-vote model 
on the dual heptagonal lattice. As above, we first present 
our numerical results for magnetization, susceptibility 
and reduced fourth-order cumulant as a function of the 
noise parameter (Fig. 0|. The qualitative properties of 
these quantities as functions of q are similar to the obser- 
vations in Figs. 0] and 01 Also here, from the intersection 
of the curves in Fig. 0Kc), we obtain q c — 0.087(5). The 
critical exponents 1/v, (3/v, and "f/v are estimated in 
Figs.Eta), (b) and (d) to 0.32(9), 0.093(6), and 0.761(9) 
respectively. The data collapse of the magnetization in 
Fig. EKc) confirms these measurements. Also here these 
values do not give a good data collapse of x m the region 
Q < Vc- 

Up to now, we obtained the critical exponents 1/v, 
/3/v and 7/f for the majority- vote model on the heptag- 
onal and dual heptagonal lattices. These exponents, and 
thus the universality classes, are different from the same 
model on square lattices. Interestingly, by checking the 
hyperscaling relation among the exponents (27| 

2(3/v + j/v = D efi , (8) 

we find that D cff = 0.125 + 0.824(7) = 0.949(7) for 
the square lattice with free boundary condition, D e g = 

0. 229 + 0.721(6) = 0.950(6) for the heptagonal lattice, 
and Doff = 0.187(2) + 0.761(9) = 0.949(1) for the dual 
heptagonal lattice. Without the finite-size effects men- 
tioned above, these values are, we believe, consistent with 

1. (If the scaling variable is chosen to be L = N 1 / 2 , 
this means D e g = 2.) In terms of these values and tak- 
ing the finite-size scaling effect into account, we propose 
that our critical exponents for the majority-vote model 
on negatively curved surface also satisfy a hyperscaling 
relation like the Rushbrooke and Josephson scaling laws 
2/5 + 7 = vd = Dcff, where d is a dimension of the un- 
derlying lattice. Since the value 1/v for both the heptag- 
onal and dual heptagonal lattices are smaller than 0.5, 
our results can be explained by an effective dimensional- 
ity of the two hyperbolic lattices greater than 2. (This 
is reasonable — the dimensionality of hyperbolic surfaces, 
embedded in an Euclidean geometry, is also larger than 
two.) 

D. Effect of boundary nodes 

From our numerical results, we know that the critical 
behaviors of the majority-vote model on the heptagonal 
and dual heptagonal lattices are different. One poten- 
tial explanation comes from by the different topological 
structure of the boundary nodes. (The boundary nodes 
of the regular heptagonal lattice interacts less with the 
inner part of the surface compared with the dual heptag- 
onal case.) 

To explore this boundary effect, one straightforward 
way is to consider a (regular or dual) heptagonal lattice 
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FIG. 4: (Color online) Majority- vote model on the heptagonal lattices. Magnetization M (a), susceptibility \ (b), and reduced 
fourth-order cumulant U (c), as a function of the noise parameter q for several values of the system size N. The critical point 
q c = 0.034(8) is estimated as the point at which the different curves for different N intercept each other. 
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FIG. 5: (Color online) (a) Log-log plot of the size dependence of the maximum values of derivatives of various thermodynamic 
quantities used to determine l/u. (b)Log-log plot of the magnetization at q = q c as a function of N. The slope of the best fit 
gives (3/u = 0.114(5). (c) Data collapse of the magnetization M shown in Fig. |4ja). The exponents used for the data collapse 
l/u = 0.3, [3/u — 0.115. (d) Log-log plot of the maximum of the susceptibility as a function of N. From this we estimate the 
critical exponent j/u — 0.721(7) as the best fit of data points, (e) Data collapse of the susceptibility shown in Fig. |4fc). The 
fitting exponents /3/u — 0.722, l/u = 0.3. 
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FIG. 6: (Color online) Majority- vote model on the dual heptagonal lattices. Magnetization M (a), susceptibility \ (b), and 
reduced fourth-order cumulant U (c), as a function of the noise parameter q for several values of the system size N. The critical 
point q c — 0.087(5) is estimated as the point where the curves for different jV intercept. 
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FIG. 7: (Color online) (a) Log-log plot of the size dependence of the maximum values of derivatives of various thermodynamic 
quantities used to determine 1/u. (b) Log-log plot of the magnetization at q = q c as a function of N. The slope of the best fit 
of the points gives f3/v — 0.093(6). (c) Data collapse of the magnetization M shown in Fig. |6ja). The exponents used for the 
data collapse 1/u = 1/3, f3/u — 0.094. (d) Log-log plot of the maximum of the susceptibility as a function of N. From this 
plot we estimate the critical exponent y/u = 0.761(9). (e) Data collapse of the susceptibility shown in Fig. [6jc). The fitting 
exponents 1/u = 1/3, f3/u — 0.762. 



with a large size, and investigate the magnetization in 
the center and compare that with a heptagonal lattice 
of a smaller size. In particular, we perform simulations 
on the heptagonal lattice at level 7 (N — 5887), and 
track the magnetization in the inner part within level 
3, 4, 5, and 6 separately. Then we compare the results 
with the magnetization of Fig. Ufa) . The same procedure 
is also done for dual heptagonal lattice at level 8 (N = 
4264), but within level 4, 5, 6, and 7 respectively. The 
numerical results are summarized in Figs. [5] and [5] We 
note that the nonvanishing boundary nodes play different 
roles in the ordering processes. As can be observed in 
Fig. [8j the magnetization in all the inner parts of a big 
heptagonal lattice are greater than that on a heptagonal 
lattice with same size, which indicates that the existence 
of a higher level of boundary can facilitate the ordering 
of the inner spins. This picture is changed, however, 
on the dual heptagonal lattice. In the subcritical region 
q < q c , the boundary nodes promote the ordering process 
as before. In the supercritical region q > q c , on the 
contrary, they impose opposite influence on the ordering 
process driving the system towards more disorder [Fig. [5]. 
A structural cause for this phenomenon is, we believe, the 
different local topologies of the two lattices. Looking at 
Fig. [TJ the nodes one step from the boundary have fewer 
connections to the boundary in the heptagonal lattice 
than its dual (the path-length going from the periphery to 
the center, for systems of the same size, are longer). The 
boundary is also more indirectly coupled to the interior in 
the heptagonal lattice. The boundary nodes of the dual 
heptagonal lattice should thus have a stronger influence 
on the configuration of the interior (we hesitate to say 



"bulk properties" since the boundary is a finite fraction 
of the interior). We expect similar stronger boundary 
effects for dual heptagonal lattices also exists for other 
statistical spin models defined. 
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FIG. 8: (Color online) Magnetization M as a function of the 
noise parameter q. (a) Filled squares represent the magne- 
tization at the inner level 3 of the heptagonal lattice with 7 
levels in total. Open circles represent the heptagonal lattice 
with 3 levels, (b), (c) and (d) show the same situation as (a) 
but for different levels. 
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mentary support for the results obtained in Ref . [20| . In 
another recent work [l|, Shima et al. studied the Ising 
model on heptagonal lattices, and found the critical expo- 
nent j/v = 0.655. Since our critical exponent 7/1/, also 
considering the estimated error, is quite different from 
this value, we believe that the majority-vote model on 
the heptagonal lattice belongs to a different universality 
class as the Ising model on this topology. This result 
can be regarded as another evidence for the conjecture 
by Yang et al. 20] mentioned above. 



IV. CONCLUSIONS 



FIG. 9: (Color online) Magnetization M as a function of the 
noise parameter q. (a) Filled squares symbolize the magneti- 
zation in the inner level 4 of the dual heptagonal lattice with 
8 levels in total. Open circles represent the heptagonal lattice 
with 4 levels, (b) , (c) and (d) show plots corresponding to (a) 
for other levels. 



E. Comparisons with the Ising model 

In Ref. [l4[ , Oliveira found that the critical exponents 
for the majority-vote model on square lattice with pe- 
riodic boundary condition are v — 0.99 ± 0.05, 7/1/ = 
1.73 ± 0.05, and (ijv = 0.125 ± 0.005. These results 
demonstrate the majority-vote model on a torus has 
the same universal critical behavior as the equilibrium 
Ising model. Recently, Yang et al. investigated a slightly 
different majority-vote model on d-dimensional hyper- 
cubic lattices. Their model is identical to the origi- 
nal one by simple replacement of the temperature pa- 
rameter to a noise parameter like in the present paper, 
tanh(l/fcT) = (1 - 2q) [2(j. Their simulation results 
suggest that the critical exponents for the majority-vote 
model in three dimensions are different from those of the 
Ising model, and that the results for four and five di- 
mensions are far from the standard mean- field values. 
For their version of the majority-vote model in two di- 
mensions, the global and local configuration energy dif- 
ferences between before and after spin-flipping are not 
identical [3(3] , but the sign of both energies is the same, 
whereas for d > 3, the sign of the two types of energy 
differences is not always the same. For the Ising model, 
however, the differences of the global and local configura- 
tion energy are exactly identical, regardless of dimension. 
They conjectured that the discordance of the sign of the 
energy difference between the global and local energy is 
responsible for the different critical behaviors. Observe 
that our critical exponents \Jv « 0.3 and 0.33 in Fig. 02a) 
and Fig. C0a), which are smaller than 0.5, so from the 
hyperscaling relation vd=\ (d is the dimension) we see 
that our hyperbolic lattices can be described as having 
an effective dimensionality greater than 2, but smaller 
than 4. In this sense, our present work provide comple- 



To summarize, we have studied the critical behavior of 
the majority-vote model on the heptagonal and dual hep- 
tagonal lattices. These lattices possess a peculiar prop- 
erty: the ratio of the size of the boundary to the to- 
tal size remains finite even in the thermodynamic limit. 
Finite-size scaling analysis reveals that the critical expo- 
nents for magnetization and susceptibility deviate from 
those of the majority-vote model on a torus (which be- 
longs to the same universality class as the equilibrium 
Ising model) and are also different from those of the 
Ising model on heptagonal and dual heptagonal lattices. 
In particular, the best fit of these exponents provided 
Pjv = 0.114(5), 7/1/ = 0.721(7) for the heptagonal lat- 
tice, and Pjv = 0.093(6), 7/1/ = 0.761(9) for the dual 
heptagonal lattice. By comparing to the majority-vote 
model on square lattices with free boundary conditions, 
we found that the free boundaries result in strong finite- 
size scaling effect, which in turn leads to the measured 
effective dimensionality smaller than unity. Nonetheless, 
we believe that the critical exponents of the majority- 
vote model defined on negatively curved surface also sat- 
isfy the hyperscaling relation ipjv + 7/1/ = D e g. 

Furthermore, for the hyperbolic lattices, we also inves- 
tigated the effect of the boundary nodes on the ordering 
process. It was shown that the boundary nodes have dif- 
ferent functions in our two lattices. For the heptagonal 
lattice, the boundary has a positive influence on order- 
ing, whereas for the dual heptagonal lattice the boundary 
nodes can either facilitate or inhibit ordering depending 
on the magnitude of the noise. These findings are further 
evidence that the underlying geometric structure deter- 
mine the critical properties of the majority-vote model. 
For the future, it will be interesting to explore the effect 
of boundary nodes on the ordering processes of other sta- 
tistical spin models defined on hyperbolic lattices. 
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